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Within recently proposed approach aimed to effectively describe the observed non-Bose type 
behavior of the intercept A of two-particle correlation function C(p,K ) of identical pions or kaons 
detected in heavy-ion collisions, the g-deformed oscillators and g-Bose gas picture are employed. 
For the intercept A, connected with deformation parameter q, the model predicts a fully specified 
dependence of A on pair mean momentum K. The intercepts A^- and A k for pions and kaons, differing 
noticeably at small K, should merge at K large enough, i.e., in the range |Kj >800 MeV/c, where 
the effect of resonance decays is negligible. In this paper we confront, fixing q appropriately, the 
predicted dependence = A^-(K) with recent results from STAR/RHIC for 7r _ 7r _ and 7r + 7r + pairs, 
and find nice agreement. Using the same q, we also predict behavior of A for kaons. 

PACS numbers: PACS numbers: 25.75.-q, 25.75.Gz, 03.75.-b 


Two-particle correlations in momentum space can be 
used to extract information about the space-time struc¬ 
ture of the emitting sources created in heavy ion col¬ 
lisions. The method exploits in an essential way the 
quantum mechanical uncertainty relation between coor¬ 
dinates and momenta, and thus any formal treatment of 
two-particle correlations must be based on a quantum 
mechanical description. For so-called ’’chaotic” sources 
where the two particles are emitted independently the 
description can be based on the single-particle Wigner 
density S(x,K) of the source (source function). 

In standard quantum mechanical treatment the Bose- 
Einstein correlations are due to symmetrization of the 
two-particle (many-particle) wave function (suppose par¬ 
ticles are emitted independently) ^> 7a7b (x a , x;,, t) = 
^['0 7a (x a ,t)^ 76 (x b ,t) + e la if la (x b , t) ip lb (x a , f)] with 
a = 0 (a = 7r) for identical bosons (fermions). The 
indices 7 a ,7f) of the 1-particle wave functions label com¬ 
plete sets of 1-particle quantum numbers. In the follow¬ 
ing, we consider two-particle correlations of noninteract¬ 
ing spin zero identical bosons. The correlation function, 
with Pi (k) and P 2 (k 0 , k b ) being single- and two-particle 
probabilities to detect particles with given momenta, is 
defined as 


G(k a ,k b ) 


P 2 (k a ,k b ) 
Pl(k a ) Pi (kb)' 


( 1 ) 


In the absence of final state interactions (FSI, see Q), for 
chaotic source, the correlation function can be expressed 
as follows ||: 


C(k a , k b ) = 1 


cos a 


I fd 4 


D S(x,K ) 


/ d A xS(x,k a ) J d 4 yS(y,k b ) 


( 2 ) 


with the 4-momenta K = f(k a + kb) as pair mean mo¬ 
mentum and p = k a — kb as relative momentum. The 


source function S(x,K) is defined by the single-particle 
states , 0 7 (x) at freeze-out time and the source density 
matrix p 77 / as, e.g., in j|j. Obviously, from (|j) at 
zero relative momentum k a = k b one gets C(k a ,k a ) = 
1 + cos a = 1 + A and since for bosons a = 0 it fol¬ 
lows that (7(k 0 , k 0 ) = 2, i.e., A = 1. To fit experimental 
data, the correlation function of identical bosons is usu¬ 
ally presented as C{p,K) = 1 + A f(p,K), with f{p,K) 
commonly taken as Gaussian so that /(p = 0, A") = 1. 
From the very first experiments it was deduced that A 
is lesser than one, the typical experimental values being 
A = 0.4 - 0.9 . The second term in (|^) is obviously due 
to quantum-mechanical interference, and deviation of A 
from unity manifests weakening of the interference ef¬ 
fects which can occur due to different reasons - influence 
of long lived resonances, coherent emission, etc. 

Let us explain the key idea of the model developed in 
|^] (named AGI-model in what follows) and further 
exploited in this letter. In two-boson correlations the 
deviation of intercept A from unity, besides the contribu¬ 
tion due to effects from long-lived resonances, can also be 
caused by the averaged softening of quantum-statistical 
effects in the peculiar short-lived many-particle systems 
formed in relativistic heavy ion collisions. In such small 
system, the symmetrization angle a of ip-y a -y b (x Q , x b , t) 
can be distorted by an additional phase due to a nonho- 
mogenuity of the system at freeze-out times (strong radial 
and azimuthal flows). These peculiarities can cause the 
effect analogous to Aharonov-Bohm one. As result, a fi¬ 
nite value of averaged symmetrization angle may appear: 
a > 0 for bosons and a < ir for fermions. 

Now, trying to explain experimental data with formula 
(B]), it is natural to relate the parameter A with the av¬ 
eraged angle a to get the reduction factor A by means of 
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cos cf. That is, the deviation of intercept A from unity is 
viewed to be due to fluctuations of symmetrization angle 
a , i.e. 


A = cos a. (3) 

Notice that slow bosons (pions, kaons) will experience 
bigger fluctuations (deviations) of symmetrization angle 
a than the particles with high velocities in fireball frame. 
That is, the deviation of intercept A from unity for slow 
bosons should be more sizable than for the fast ones. 

To implement our key idea we exploit quantum field 
theory with g-deformed commutation relations (qDCR) 
and the techniques of g-boson statistics (see [Q and refs, 
therein) which reflects a partial suppression of the quan¬ 
tum statistical effects. In || || it was argued that the 
algebra of qDCR is connected, for real q only, with 
the so-called nonextensive statistics introduced by Tsal- 
lis JtJ. This type of generalized statistics has already 
found numerous applications in diverse branches of mod¬ 
ern physics (see [|j for refs.). In particular, nonexten¬ 
sive statistics was applied to the problems of high-energy 
nuclear collisions (||] and refs, therein). However, the 
techniques of g-boson statistics based on qDCR allows 
the use of complex values as well as the real values for 
the deformation parameter q , depending on the choice 
of algebraic realization of qDCR. The physical reasons 
for usage of qDCR, and subsequent interpretation of q, 
essentially differ depending on the case of q real or com¬ 
plex. Introducing deformed statistics with q real enables 
one to effectively account for interaction effects by means 
of non-interacting ideal gas of “modified” particles. On 
the other hand, the approach based on qDCR provides 
the ability to model the effects involving the Aharonov- 
Bohrn like phase, intimately connected with symmetriza¬ 
tion properties of wave functions. 

For the system of pions or kaons produced in heavy ion 
collisions, we employ the ideal g-Bose gas picture. Physi¬ 
cal meaning or explanation of the origin of g-deformation 
in the considered phenomenon sharply differs in the case 
of real deformation parameter g from the case when g is 
a pure phase factor, as will be seen in what follows. 

The AGI-model exploits two different sets of qDCR. 
The first is the multimode Biedenharn-Macfarlane (BM- 
type) g-oscillator defined as ju|: [. Nj,bj\ = —bj, 
[Nj,b^] = &t, bjb^ — q~ 1 b^bj = q Nj , where different 
modes (i ^ j) commute. Then, b\bi = [iV,] g (here the 
“g-bracket” means [r] q = ( q r — q~ r )/{q — g _1 ) ) so that 
b\bi = Ni is recovered in the “classical” (“no deforma¬ 
tion”) limit g —> 1 . Below, for the BM-type g-oscillators 
it is meant that 

g = exp(?'6 l ) , 0 < 9 < 7r/2 . (4) 

The second multimode g-oscillator used in AGI-model 
is the set of Arik-Coon (AC-type) g-oscillators, defined 


by the relations [JTlJ [N, a] = —a, [N, a^] = a*, and 
aa t — g<Aa = 1 (subscript suppressed). Again, at g^l, 
the bilinear aj ai does not equal the number operator A/) 
(as is true for usual bosonic oscillators, i.e., at g = 1). 
Instead, afa^ = [[A/)]] where now the notation [[r]] = 
(1 — g r )(l — g), is used. The g-bracket [[A]] for an oper¬ 
ator A is understood as a formal series. At g —> 1, from 
[[A]] one recovers A. In what follows we set 

- 1 < g < 1 . (5) 


For each such value of the deformation parameter g, the 
a\ and ai are mutual conjugates. Note that the in¬ 
verse of the relation a\ai = [[A/)]] is given by a for¬ 
mula expressing the operator A/) as a formal series of 
creation/annihilation operators. 

For a multi-pion (-kaon) system, viewed as ideal gas of 
g-bosons the Hamiltonian is taken as 


H = J2“iN x ( 6 ) 

i 

with i labelling energy eigenvalues, w* = m 2 + k.f , and 
Ni defined as above. This is unique truly noninteracting 
Hamiltonian with additive spectrum [|J . We assume dis¬ 
crete 3-momenta of particles (the system is in a box of 
volume ~ L 3 ). For the set of AC-type g-oscillators, one 
takes Ni instead of N t in (](]). 

Statistical properties are obtained by evaluating ther¬ 
mal averages (A) = Sp(Ap)/Sp(p) , p = e~^ H , with the 
Hamiltonian (y) and (3 = 1/T. 

With b\bi = [Ni] q and q + g” 1 = [2] g = 2cos0, the 
g-deformed distribution function is obtained as I! lil] 


(b\bi) 


1 

e /3wi _ 1 + Si ’ 


S,=2 


1 — cos 8 
1 _ e -/3uJt ■ 


( 7 ) 


If 9 —> 0, it yields Bose-Einstein (BE) distribution. Note 
that the g-distribution function (0) is real. 

The g-distribution (f\) deviates from the quantum 
Bose-Einstein just in the “right direction” towards the 
classical Boltzmann distribution, that reflects a decreas¬ 
ing of quantum statistical effects. For kaons, whose mass 
itik is bigger than m„, analogous curve should lie closer, 
than pion’s one, to that of BE distribution [||. 

In the case of AC-type g-bosons with real g from (£j), 
one arrives at the distribution function (cf. 

<°i a i) = e 0^q ■ 


In the no-deformation limit g —> 1, this also reduces to 
the Bose-Einstein distribution, since at g = 1 we return 
to the standard system of bosonic commutation relations. 

The deviation from standard BE statistics is natural 
thing if one considers the system of interacting particles 
versus that of non-interacting particles (ideal gas). For 
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instance, natural type of interaction is the hard-core re¬ 
pulsion of the particles that assumes particle finite self¬ 
volume. This type of interaction, as was shown in [ p~3| , 
results in the same kind (|t|) of modified statistics. At mi¬ 
croscopical level, a finite self-volume arising due to com¬ 
posite structure of particles results in g-deformed commu¬ 
tation relations Q and subsequently results in certain 
g-deformed statistics of the gas of such particles. 

The two-particle distribution corresponding to the 
BM-type g-oscillators is 


(b\b\bA) 


2 cosf? 

e 2(3ui _ 2 cos(2 6)eP Ui + 1 


(9) 


From this and eq. (§, one obtains the intercept A = 
Aj + 1 = (b\b\bibi)/((b\bi)) 2 of two-particle correlations 
(subscript omitted) as 


2cos0 (cosh(/3w) — cost?) 2 
(cosh (Plo) — 2 cos 2 8 + l)(cosh(/3w) — 1) ’ 


At (3u) —> oo (i.e., at low temperature and fixed momenta 
or large momenta and fixed temperature) the asymp¬ 
totics of intercept is given merely by the deformation 
angle 8 (recall that q = exp(i0)): 


A = A asymp = 2cos6»-l (T —> 0 or |K| -> oo). (11) 

From this and eq. (|J) we have the (asymptotical) relation 
cos 8 = cos 2 Note that, if the unique cause forcing the 
intercept to be lesser than one is the decays of resonances 
(the conventional viewpoint), all the curves would tend 
to the value A = 1 in the large |K| limit. In contrast, we 
predict a constant A < 1, as in (|ll|). 

In the case of AC-type g-oscillators, the formula 
(a\a\aiai) = (1 + q){e ISui — g) _1 (e /3a ’ i — g 2 ) -1 for two- 
particle distribution combined with (j?]) leads to 

A = — 1 + (a) a) aa) / (a) a) 2 = q — ^-1-(12) 

e"“ — q z 

In this case, for T —> 0 or K| —> oo we have A asymp = g. 

Below, the two versions ([bjj) and corresponding 
to BM- and AC-types of g-deformation are compared to 
the recent STAR/RHIC data. The experimental values 
for intercept parameter A in Figs. [l| and ^ are taken from 
Ref. [^5| . The theoretical values are obtained by averag¬ 
ing over given rapidity y and transverse momentum K t 
intervals Aj = K{’ max — Jv/’ min , j = 1,2,3: 


A y /2 K{ 

A j = -£- J J dK t Hq,m,T,y,K t ), (13) 

V -A v /2 ^ K i,min 

where m is particle mass. 

Expressions (0) and & for A(g, to, T, y, K t ) were 
used in (0) for obtaining theoretical points shown in 
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FIG. 1: The transverse momentum |K t | dependence of the 
intercept A of two-pion correlation. Deformation parameter q 
is a real quantity, 0 < q < 1. 
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FIG. 2: The transverse momentum |Kt| dependence of the in¬ 
tercept A of two-pion correlation. The deformation parameter 
q is taken in the form q = e t8 . 


Fig. 0 and Fig. |] respectively. One can see from these 
figures that the agreement of experimentally measured 
values of intercept parameter A with the theoretically 
calculated ones is very good. 

Detailed comparison with the experiment []l5| gives: 
the values A,; obtained from (13) at real g, see (12), fit 
better the three experimental values for the intercept of 
7r + 7r + correlations. On the other hand, the values calcu¬ 
lated by (13) with g a pure phase factor, see (10), agree 
better with the three experimental values for the inter¬ 
cept of 7r - 7r“ correlations. Possible explanation of the 
observed difference between experimental values of inter¬ 
cept for 7r _ 7r _ -pairs and 7r + 7r + -pairs could be the influ¬ 
ence of the Coulomb FSI of these charged pions with the 
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FIG. 3: The transverse momentum |Kt| dependence of the 
intercept A of two-kaon correlation. Both the case of real 
q = 0.63 (lower triple of curves) and the case of q = e lB 
(upper triple of curves) with 9 = 28.5° are shown. 


positive charge of fireball protons. Since the AGI-model 
predicts that parameter A will asymptotically reach a 
constant value A asymp ’ < 1, determined by q only, at 
sufficiently large (500 — 600 MeV/c) pion pair mean mo¬ 
mentum | K |, in order to check this prediction measure¬ 
ments at higher Kt are necessary. Such measurements 
should be available in near future at RHIC. 

For the prediction of intercept of kaons, we will use 
the values of q which provide the best fit of experimental 
data for pions (see Figs. |j] and ||): q = 0.63 or 9 = 28.5°, 
assuming a universality of the deformation parameter for 
description of excited hot hadronic matter. The result of 
averaging in rapidity —0.5 < y < 0.5, given by first in¬ 
tegral in ©> is shown in Fig. as solid curve in each 
of the triples of curves. The other two curves in each 
triple correspond to fixed value of rapidity: y = 0 (dot¬ 
ted curve) and y = 0.5 (dashed curve). Note that y = 0 
curve and solid curve almost coincide. As it is clearly 
seen, the cases of real q and q a phase factor supply sig¬ 
nificantly different values for the kaon intercept A k- It is 
tempting to use just this feature for making preference of 
particular version of the deformation parameter q - real 
or pure phase. The choice is important because differ¬ 
ing physics is behind these two versions: real q may re¬ 
flect for instance particle finite size effects |li| or particle 
composite structure [ pdf , and complex valued q may refer 
to deformed symmetrization properties of wave functions 
(like in Aharonov-Bohm effect), relevant for short-lived 
systems occurred in heavy-ion collisions. It is also pos¬ 
sible that a phase-type q encodes [jib! the effects from 
mixing at the composite (quark) level. Recent data from 
NA44 H|, i.e. A = 0.84 ±0.13 and A = 0.61 ±0.36 resp. 
for (K t ) ss 0.25 GeV/c and (K t ) « 0.91 GeV/c do not 
yet help in making choice of optimal version for q. 

In summary, we have presented comparison of the AGI 


model with experimental data on two-particle correla¬ 
tions at RHIC, and found remarkable agreement. We 
used the parameters extracted from comparison with 
pion’s data to predict behavior of intercept of kaon corre¬ 
lation functions. We stress again the crucial importance 
of correlation measurements at high transverse momenta 
in order to check the predicted asymptotical ’'saturation'’ 
of intercept parameters. Measurements |K| in the range 
up to 500 — 600 MeV/c for pions (up to 700 — 800 MeV/c 
for kaons) should be possible by RHIC detectors such as 
STAR and PHENIX. The asymptotical behavior of in¬ 
tercept parameter A within the proposed model, see 0 ) 
for phase-type q, should determine the actual value of 
the deformation parameter q supposed to be a universal 
quantity for relativistic heavy ion collisions. 
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